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Abstract. We study the effect of screening by bound electrons in low energy nuclear reactions. We use 
molecular dynamics to simulate the reactions involving many electrons: D+d, D+D, ^He+d, '^He+D, ^Li+d, 
^Li+D, '^Li+p, '^Li+H. Quantum effects corresponding to the Pauli and Heisenberg principles are enforced 
by constraints in terms of the phase space occupancy. In addition to the well known adiabatic and sudden 
limits, we propose a new "dissipative limit" which is expected to be important not only at high energies 
but in the extremely low energy region. The dissipative limit is associated with the chaotic behavior of 
the electronic motion. It affects also the magnitude of the enhancement factor. We discuss also numerical 
experiments using polarized targets. The derived enhancement factors in our simulation are in agreement 
with those extracted within the i?-matrix approach. 



PACS. 25.45.-z - 34.10.+X 

1 Introduction 

The relation between the tunneling process and dynam- 
ical chaos has been discussed with great interests in re- 
cent years. PI2] Though the tunneling is completely quan- 
tum mechanical phenomenon, it is influenced by classical 
chaos. In the sense that the the chaos causes the fluctua- 
tion of the classical action which essentially determines the 
tunneling probability. We study the phenomenon by ex- 
amining the screening effect by bound electrons in the low 
energy fusion reaction. In the low energy region the exper- 
imental cross sections with gas targets show an increasing 
enhancement with decreasing bombarding energy with re- 
spect to the values obtained by extrapolating from the 
data at high energies 0. Many studies attempted to at- 
tribute the enhancement of the reaction rate to the screen- 
ing effects by bound target electrons. In this context one 
often estimates the screening potential as a constant de- 
crease of the barrier height in the tunneling region through 
a fit to the data. A puzzle has been that the screening po- 
tential obtained by this procedure exceeds the value of the 
so called adiabatic limit, which is given by the difference of 
the binding energies of the united atoms and of the target 
atom and it is theoretically thought to provide the maxi- 
mum screening potential 0]. Over these several years, the 
redetermination of the bare cross sections has been pro- 
posed theoretically and experimentally 6 , using the 
Trojan Horse Method |7IIHll21| . The comparison between 
newly obtained bare cross sections, i.e., astrophysical S- 
factors, and the cross sections by the direct measurements 



gives a variety of values for the screening potential. These 
values are often smaller than the sudden limit or larger 
than the adiabatic limit. Theoretical studies performed us- 
ing the time-dependent Hartree-Fock(TDHF) scheme 
[TU] suggest that the screening potential is between the 
sudden and the adiabatic limits. 

One of the aims of this paper is to try to assess the 
effect of the screening quantitatively. Up to now, the dy- 
namical effects of bound electrons have been studied only 
in some limited cases with a few bound electrons(the D+d 
with atomic target pilIU| and molecular D2 target JT] , the 
•^He-l-ii |9j) with the TDHF method. We investigate here 
the dynamical effects, including the tunneling region, for 
other systems with many bound electrons; D-l-D, '^He-I-D, 
looking the effect of the electron capture of projectile. 
We see also some reactions including Li isotopes; ^Li+d, 
^Li+D, '^Li+p and '^Li+H. 

To simulate the effects of many electrons, we use the 
constrained molecular dynamics (CoMD) model pill2lll5j . 
At very low energies fluctuations are anticipated to play a 
substantial role. Such fluctuations are beyond the TDHF 
scheme. Not only TDHF calculations are, by construc- 
tion, cylindrically symmetric around the beam axis. Such 
a limitation is not necessarily true in nature and the mean 
field dynamics could be not correct especially in pres- 
ence of large fluctuations. Molecular dynamics contains 
all possible correlations and fluctuations due to the initial 
conditions (events). For the purpose of treating quantum- 
mechanical systems like target atoms and molecules, we 
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use classical equations of motion with constraints to sat- 
isfy the Heisenberg uncertainty principle and the Pauli 
exclusion principle for each event 12 . In extending the 
study to the lower incident energies, we would like to stress 
the connection between the motion of bound electrons and 
chaos. In fact, depending on the dynamics, the behavior of 
the electron(s) is unstable and influences the relative mo- 
tion of the projectile and the target. The feature is caused 
by the nonintegrablility of the A^-body system(Af > 3) 
and it is well known that the tunneling probability can 
be modified by the existence of chaotic environment. We 
discuss the enhancement factor of the laboratory cross sec- 
tion in connection with the integrability of the system by 
looking the inter-nuclear and electronic oscillational mo- 
tion. More specifically we analyze the frequency shift of 
the target electron due to the projectile and the small os- 
cillational motion induced by the electron to the relative 
motion between the target and the projectile. We show 
that the increase of chaoticity in the electron motion de- 
creases the fusion probability. 

The paper is organized as follows. In sect. El we intro- 
duce the enhancement factor /e and describe the essence of 
the Constrained molecular dynamics approach briefly. In 
sect.Elwe apply it to asses the effect of the bound electrons 
during the nuclear reactions. We discuss also the relation 
between the amplitudes of the inter-nuclear oscillational 
motion and the enhancement factor. We summarize the 
paper in sect.^ 



2 Formalism 

2.1 Enhancement Factor 

We denote the reaction cross section at incident energy in 
the center of mass E by (t{E) and the cross section ob- 
tained in absence of electrons by ctq (E) ■ The enhancement 
factor /e is defined as 



fe 



ME)' 



(1) 



If the effect of the electrons is well represented by the 
constant shift Ue of the potential barrier, following |13ll^ . 
{Ue < E): 



fe ~ exp 



E 



(2) 



where ri{E) is the Sommerfeld parameter |14j . 



2.2 Constrained Molecular Dynamics 

We estimate the enhancement factor fe numerically using 
molecular dynamics approach; 



dvi PiC^ dp 



dt £i 



(3) 



where (ri,pi) are the position, momentum of the parti- 
cle i at time t. Ei = ^/pfc^ + rn 



C/(r,;) and mi are 



its energy. Coulomb potential and mass, respectively. We 
set the starting point of the reaction at lOA inter-nuclear 
separation. In Eqs. Q we do not take into account the 
quantum effect of Pauli exclusion principle and Heisen- 
berg principle. As it is well known that these classical 
equations lO can be derived by using the variational cal- 
culus of Lagrangian C of the classical system as well. So 
as to take the feature of the Pauli blocking into account 
in this framework, we use the Lagrange multiplier method 
for constraints. 

Our constraints which correspond to the Pauli block- 
ing is /i < 1 in terms of phase space density, note that the 
phase space density can be directly related to the distance 
of two particles, i.e., r^py , in the phase space. Here Vij = 
\ri — rj \ and — \pi~pj\. The relation /i < 1 is fulfilled, 
if TijPij > £,ph6si,Sj, where = 2tt{3/4:Tt)'^/^. i,j refer 
only to electrons and Si,Sj{= ±1/2) are their spin pro- 
jection. For the Heisenberg principle VijPij > S,h^, where 
f// = 1, j and j refer to not only electrons but the nucleus. 
It is determined to reproduce the correct energy of hydro- 
genic atoms. Obviously the conditions r^py = £,h{p)^ 
must be fulfilled in the ground state configuration rather 
than TijPij > £,H{p)'h. 

Using these constraints, the Lagrangian of the system 
can be written down as 



E ( ^^5.5, - 1 ) , (4) 



where Af and are Lagrange multipliers. The varia- 
tional calculus leads 



dt 



Pi(_ 



-y 



' dpi ' 



(5) 



4h 4p 



9r. 



dr. 



In order to obtain the atomic ground-state configuration. 
We perform the time integration of the eqs. lO and JSJ. 
The value of and Af are determined depending on 
the magnitude of VijPij. If r^py is (smaller)larger than 
^H{P)^j ^ h^s positive(negative) sign. Thus we change the 
phase space occupancy of the system. The constraints re- 
strict us to variations AC = that keep the constraints 
always true |15|. In this way we obtain many initial con- 
ditions which occupy different points in the phase space 
microscopically. 

In order to treat the tunneling process, we define the 
collective coordinates R'^o" and the collective momentum 

pco// 



as 



■\coll 



Pp Pt, 



(7) 



where rT,rp (pt,Pp) are the coordinates (momenta) of 
the target and the projectile nuclei, respectively. When 
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the collective momentum becomes zero, we switch on the 
collective force, which is determined by Fp " = P'=°" and 
F™" = — pt^o''^ to enter into imaginary time ^Hl- We fol- 
low the time evolution in the tunneling region using the 
equations, 



dr. 



't(p) _ Pt(p) 



dp 



T(P) 



dr 



& 



T{P) 



dr 



VrC/(r?(p))-2F 



T(P)' 



, (8) 

where r is used for imaginary time to be distinguished 
from real time t. r^(p) and Y*t(p) ^'^^ position and mo- 
mentum of the target (the projectile) during the tunnel- 
ing process respectively. Adding the collective force corre- 
sponds to inverting the potential barrier which becomes 
attractive in the imaginary times. The penetrability of the 
barrier is given by |16| 



n{E) = {l + ciip{2AiE)/h)y 
where the action integral A{E) is 

A{E) ^ I P'^"" dR""" 



(9) 



(10) 



and rh are the classical turning points. The internal 
classical turning point r?, is determined using the sum of 
the radii of the target and projectile nuclei. Similarly from 
the simulation without electron, we obtain the penetrabil- 
ity of the bare Coulomb barrier IIq{E). 

Since nuclear reaction occurs with small impact pa- 
rameters on the atomic scale, we consider only head on 
collisions. The enhancement factor is thus given by eq. (|T|l. 



/e = n{E)/no{E) 



(11) 



for each event in our simulation. Thus we have an ensemble 
of fe values at each incident energy. 



3 Application to the Electron Screening 
Problem 

3.0.1 D+d and D+D reactions 

Fig. ^ shows the incident energy dependence of the en- 
hancement factor for the reactions D-l-d and D-l-D, where 
the systems involve 1 and 2 electrons respectively. The 
open and closed squares show the average enhancement 
factors fe over events for the reactions D+d and D+D, 



respectively. The variances U = y /| — (/e)^ are shown 
with error bars. The dotted and dash-dotted curves show 
the enhancement factors in the adiabatic limit /i"^^'' for 
an atomic deuterium target and it is obtained by assuming 
equally weighted linear combination of the lowest-energy 
gerade and ungerade wave function for the electron, re- 
flecting the symmetry in the D+d, i.e.. 



(AD) 



1 I f.^r,(E)- 



+ e 



(12) 



CD 



I I 1 — I 1 — I 

D+d I — B — I 
I A D+D 
AD(D+d) 

DL(D+d) 

AD(D+D) 

DL(D+D) 
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Fig. 1. Enhancement factor as a function of incident center- 
of-mass energy for the reactions D-l-d and D+D. Error-bars 
represent the variances obtained from the events generated for 
each beam energy. 



where t/i^^ = 40.7 eV and C/i"^ = 0.0 eV [TUIIU] for D+d 
case. If we take into account the electron capture of the 
projectile, i.e., in the case of D+D, the enhancement factor 
in the adiabatic limit is 



(1") 



(13) 



where C/i^ "'^ = 51.7 eV and Uli"''' = 31.9 eV [HI. The 
solid curve and dashed curve show the enhancement fac- 
tors in the dissipative limit fi^^^ for the reactions D+d 
and D+D respectively. Notice how the calculated enhance- 
ment factor with their variances nicely ends up between 
the adiabatic and the dissipative limits. We performed also 
a fit of our data using eq. ^ including the very low energy 
region and obtained Ue = 15.9 ± 2.0 eV for D+d case and 
Ue =21.6 ± 0.3 eV for D+D. 

Now we look at the oscillational motions of the parti- 
cle's coordinates as the projection on the z-axis (the re- 
action axis). We denote the ^-component of rT,rp and 
Te as ZT,zp and Ze, respectively. Practically, we examine 
the oscillational motion of the electron around the target 
ZTe = Ze — zt and the oscillational motion of the inter- 
nuclear motion, i.e., the motion between the target and 
the projectile, Zs = zt + zp, which essentially would be 
zero due to the symmetry of the system in the absence 
of the perturbation. In Fig. |21 these two values are shown 
for 2 events, which have the enhancement factor fe — 
170.8 (ev. A), and fe = 6.5 (ev. B), at the incident energy 
Ecm — 0.15 keV. The panels show the Zs, ZTe as a function 
of time. The stars indicate the time at which the system 
reaches the classical turning point. It is clear that in the 
case of event B the orbit of the electron is much distorted 
from the unperturbed one than in event A. Characteris- 
tics of Zs are that (1) its value often becomes zero, as it is 
expected in the un-perturbed system, and (2) the compo- 
nent of the deviation from zero shows periodical behavior. 
It is remarkable that the amplitude of the deviation be- 
comes quite large at some points in the case of event B 
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Fig. 2. The osciUational motion of the electron around the target (lower panels) and the inter-nuclear motion (upper panels) 
as a function of time, in atomic unit, for two events, with large /<;(ev. A) and small /e(ev. B), for the D-l-d reaction at the 
incident energy O.lSkeV. The inter-nuclear separation is lOA at f = 0. 



which shows the small enhancement factor. Note that in 
event B one observes clear beats, i.e., resonances. Thus 
for two events, with the same macroscopic initial condi- 
tions, we have a completely different outcome, which is 
a definite proof of chaos in our 3-body system. We can 
understand these results in first approximation by con- 
sidering the motion of the ions to be much slower than 
the rapidly oscillating motion of the electrons. From 
the FigElwe can deduce the following important fact. If 
the motion of the electron is initially in the plane per- 
pendicular to the reaction axis, the enhancement factor is 
large, event A(notice |zTe| <C i.e., the Bohr radius, 
at t ^ 0). On the other hand if there is a substantial 
projection of the electron motion, as in event B(the am- 
plitude of \zTe\ ~ i?B at t ^ 0), on the reaction axis the 
enhancement factor is relatively small because of the in- 
crease of chaoticity. The fact suggests that if one performs 
experiments at very low bombarding energies with polar- 
ized targets, the enhancement factor can be controlled by 
changing the polarization. The largest enhancement would 
be gained with targets polarized perpendicularly to the 
beam axis. 
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Fig. 3. Incident energy dependence of the enhancement factor 
for Px and Py targets. 



bound electron gives the effect of hindrance to the tunnel- 
ing probability. 



In order to test this estimation, we prepared ensembles 
of target atoms which are polarized perpendicular(Pj^) 
and parallel(P||) to the beam axis, numerically. In Fig. 13 
we show the incident energy dependence of the average 
enhancement factor for the and Py targets with pluses 
and crosses, particularly in the low energy region. The en- 
hancement factors from the P_l targets are always larger 
than that from the Py targets. In contrast to the average 
enhancement from the P_l targets, which increases mono- 
tonically as the incident energy becomes smaller, the av- 
erage enhancement from the P|| targets fluctuates. It has 
also large variances at low energies. Remarkable thing is 
that with the parallel targets the enhancement factor of- 
ten becomes less than 1. It means that in this case the 



3.0.2 ^He+d and ^He+D reactions 

An excess of the screening potential was reported for the 
reactions '^He-l-d with atomic gas ^He target, and D2 -I- 
•^He with deuterium molecular gas target, for the first time 
in the reference [31. Since then various experiments have 
been performed for these reactions. The incident energy 
covers from 5 keV to 50 keV for ^He-l-ii. Though once 
the problem of the discrepancy between experimental data 
and theoretical prediction seemed to be solved by consid- 
ering the correct energy loss data ^^li recent measure- 
ments using measured energy loss data |20) report larger 
screening potentials than in the adiabatic limit for both 
reactions. 
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The electron capture by the projectile plays a minor 
role in the case of '^He+d, since electrons are more bound 
in helium targets. However in the recent measurement 
Aliotta et al. was performed using molecular and 
targets [3D]. Thus we assess the contribution from the re- 
action ^He+D, as well. 

The enhancement factor in the adiabatic limit give 
J7e=119 eV for ^Re+d and Ue=nO eV for ^He+D, respec- 
tively. These are shown in the figure^jwith the solid curve 
for ^He-|-d and with the dashed curve for '^He-I-D. The 
comparison of these two adiabatic limits implies that the 
electron capture of projectile would give a hindrance com- 
pared with the bare deuteron projectile. Meanwhile the 
latest analysis of the experimental data using i?-matrix 
two level fit suggests the screening potential Ue ~ 60 
eV(corresponding enhancement factor is shown with dot- 
ted curve). The comparison between direct measurement 
and an indirect method, the Trojan Horse method, sug- 
gests the screening potential Ue — 180±40 eV ( the cor- 
responding enhancement factor is shown with dot-dashed 
curve) [2] . The average enhancement factors fe over events 
in our simulations using the CoMD are shown with the 
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Fig. 4. Enhancement factor as a function of incident center- 
of-mass energy for the reactions "^He+d and ^He+D. 



open and closed squares for the reactions "^He+d and '^Hc+D, 
respectively. The enhancement factors of the both reac- 
tions '^He-f d and '^He+D are in agreement with the ex- 
tracted values using the i?-matrix approach within the 
variances over all the events. Notice that our calculated 
enhancement factors for the two systems display an oppo- 
site trend as compared to the adiabatic limits. The average 
enhancement factor of the reaction ^He-|-D agrees with the 
estimation of the adiabatic limit and the reaction ■^He-l-d 
is below the corresponding adiabatic limit. The paradoxi- 
cal feature comes from the fact that an electron between 
the two ions is often kicked out during the reaction pro- 
cess, i.e., the electron configuration seldom settles down 
the ^Li+ ground state in the reaction -^He+d. It is known 
as autoionization in the context of the Classical Trajec- 
tory Monte Carlo method • Instead in the case of the 
^He-|-D, the deuterium projectile brings its bound electron 



in a tight bound state around the unified nuclei of '^He and 
d, practically it ends up with a ground state configuration 
of the ^Li atom. The fits of the obtained enhancement fac- 
tors suggests the screening potentials Ue = 82.4 ± 1.9 eV 
for the ^Re+d and Ue = 102.8 ± 3.0 eV for the ^He-hD. 



3.0.3 '^U+d, 6Li+D, ^Li+j3 and ^Li+H 



The S-factors for the reactions ^Li-|-c?, ^Li+p and ''Li-fp 
were measured over the energy range 10 keV < Ecm < 500 
keV by Engstler,ei al. They used LiF sohd targets 
and deuteron projectiles as well as deuterium molecular 
gas targets and Li projectiles. 

In the case of LiF target which is a large band gap in- 
sulator, one often approximates the electronic structure of 
the target ^Li(^Li) state by the ^Li+(''Li+) with only two 
innermost electrons. Thus for all three reactions one ex- 
pects the screening potential in the adiabatic limit t/i'^^'' = 

371.8 - 198.2 ~ 174 eV. Instead if one uses the ground 
state of the ^Li(^Li) atom and of the bare deuteron target 

as the initial state, ui^^^ =186 eV which is given by 
the solid curve in Fig. • 

However one should be aware that the deuteron or hy- 
drogen projectile plausibly moves with a bound electron 
in LiF solid insulator target |25| . Under such an assump- 
tion we could estimate the screening potential Ue"^^^ = 

389.9 - 198.2 ~ 192 eV. In the case of molecular D2 or 
H2 gas targets, as well, we should consider the electron 
capture by the lithium projectile. 

The bare S-factors for the same reaction have been ex- 
tracted using an indirect method, the Trojan-Horse Method 
through the reaction ^Li(^Li,aQ;)''He The comparison 
between direct and the indirect methods gives the screen- 
ing potential Ue — 320±50 eV. The corresponding en- 
hancement factors are shown with the dash-dotted curve. 
The contrast between the direct measurement data and 
the theoretical estimation for the bare S-factor using the 
i?- matrix theory gives Ue=240 eV. It is shown with dotted 
line. The extracted Ue with the two different methods are 
larger than the adiabatic limit. 

We simulate the reactions ^Li-l-d, ^Li-fD, ^Li-t-p and 
^Li+H. In the figure |3J and O the open and closed squares 
show the enhancement factor for the reactions ^Li-|-d and 
^Li-|-D,(and ^Li+p and '^Li-I-H) respectively. 

Again the average enhancement factors of the reaction 
^Li-l-D(^Li-l-H) are larger than those of the ^Li-t-dC^Li-l-p). 
The enhancement factors of the reaction ^Li+D are in 
agreement with the extracted values using the i?-matrix 
approach within the variances over all the events. The fit 
of the obtained average enhancement factors suggests the 
screening potentials Ue = 152.0 ± 9.9 eV for ^Li-^d and 
Ue = 214.4±18.5 for ''Li+D. The screening potential for 
the reaction ^Li-|-d in our simulation does not exceed the 
adiabatic limit nor extracted values using the i?-matrix 
theory and THM, but one for ®Li+D verges on the ex- 
tracted values using the i?-matrix approach. 
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Fig. 5. same as Fig.|l]but for the reactions "Li+d and ®Li+D. 
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Fig. 6. same as Fig.0]but for the reactions ^Li+p and ^Li+H. 

4 Summary 

We discussed the effect of the screening by the electrons 
in nuclear reactions at the astrophysical energies. We per- 
formed molecular dynamics simulations with constraints 
and imaginary time for the reactions D+d, D+D, ■^He+rf, 
^He+D, ^Li+d, ^Li+D, "^Li+p, ^Li+H. For all the reac- 
tions it is shown that both the average enhancement fac- 
tors and their variances increase as the incident energy 
becomes lower. Using bare projectiles we obtained the av- 
erage screening potential smaller than the value in the 
adiabatic limit for all reactions. It is because of the exci- 
tation or emission of several bound electrons during the 
reactions. The comparison between bare and atomic pro- 
jectile cases for each reactions revealed that the electron 
capture of the projectile guides to larger enhancements. 
The derived enhancement factors in our simulation are in 
agreement with those extracted within the i?-matrix ap- 
proach including the variances over all the events. 

We report also the results of the numerical experiments 
using polarized targets for the reaction D-l-d. Using 
targets we obtained relatively large enhancements with 
small variances, instead P|| target gives large variances of 
the enhancement factors and relatively small averaged en- 



hancement factors. It is because with the P|| targets the 
force exerted from the electron to the relative motion of 
the nuclei is oscillational, in the direction of the beam 
axis, and the motion of the electron becomes often ex- 
cited or unstable. It is the case where the chaoticity of 
the electron motion affects the tunneling probability and 
at the same time the enhancement factor of the cross sec- 
tion. This suggests that if one performs experiments at 
very low bombarding energies with polarized targets, the 
enhancement factor can be controlled by changing the po- 
larization. The largest enhancement with targets polarized 
perpendicularly to the beam direction. 
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